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Abstract 

We show theoretically that polariton pairs with a high degree of polarization entanglement can 
be produced through parametric scattering. We demonstrate that it can emerge in coincidence 
experiments, even at low excitation densities where the dynamics is dominated by incoherent 
photoluminesce. Our analysis is based on a microscopic quantum statistical approach that treats 
coherent and incoherent processes on an equal footing, thus allowing for a quantitative assessment 
of the amount of entanglement under realistic experimental conditions. This result puts forward 
the robustness of pair correlations in solid-state devices, even when noise dominates one-body 
correlations. 
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3j. Polaritons on the other 



The concept of entanglement has played a crucial role in the development of quantum 
physics. It can be described as the correlation between distinct subsystems which cannot 
be reproduced by any classical theory (i.e. quantum correlation). It has gained renewed 
interest mainly because of the crucial role that such concept plays in quantum informa- 
tion/computation as a precious resource enabling to perform tasks that are either im- 
possible or very inefficient in the classical realm j^. Scalable solid-state devices will make 
use of local electronic states to store quantum correlations 
hand, as hybrid states of electronic excitations and light, are the most promising solution for 
generation and control of quantum correlations over longer range In particular, thanks 
to the Coulomb interaction acting on the electronic part of the polariton state, resonantly 
generated pump polaritons scatter into pairs of signal and idler polaritons, in a way that 
fulfills total energy and momentum conservation. The generated polariton pairs are ideally 
in an entangled state 6|, |7| . The peculiar energy- momentum dispersion of microcavity po- 
laritons has the advantage of allowing several configurations of parametric scattering, that 
can be easily selected by setting the frequency and angle of both the pump and the detected 
beams [sl, In order to address entanglement in quantum systems IQ, 11], the preferred 
experimental situations is the few-particle regime in which the emitted particles can be de- 
tected individually [3] • In a real system, environment always act as an uncontrollable and 
unavoidable continuous perturbation producing decoherence and noise. Even if polariton 
experiments are performed at temperature of few Kelvin [3] , polaritons created resonantly 
by the pump can scatter, by emission or absorption of acoustic phonons, into other states, 
acquiring random phase relations. These polaritons form an incoherent background (i.e. 
noise), responsible of pump-induced photoluminescence (PL), which competes with coher- 



ent photoemission generated by parametric scattering, as evidenced by experiments 



9|. As 



a consequence, noise represents a fundamental limitation, as it tends to lower the degree 
of non-classical correlations or even completely wash it out [l^ . Hence understanding 
the impact of noise on quantum correlations in semiconductor devices, where the electronic 
system cannot be easily isolated from its environment, is crucial. 

In this letter, we present a microscopic study of the infiuence of time-dependent noise 
on the entanglement of polaritons generated in parametric PL. Our treatment accounts for 
realistic features such as detectors noise background, detection windows, dark-counting etc., 
needed Il5| in order to seek and limit all the unwanted detrimental contributions. We show 
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how a tomographic reconstruction [l^ , based on two-times correlation functions, can provide 
a quantitative assessment of the level of entanglement produced under realistic experimental 
conditions. In particular, we give a ready-to-use realistic experimental configuration able 
to measure the Entanglement of Formation [l^], out of a dominant time-dependent noise 
background, without any need for post-processing fll |. 

Third order nonlinear optical processes in quantum well excitons (with spin a = ±1) 
can be described in terms of two distinct scattering channels: one involving only excitons 
(polaritons) with the same circular polarization (co-circular channel); and the other (counter- 
circular channel) due to the presence of both bound biexciton states and four-particle scat- 
tering states of zero ang ular momentum (J = 0) [l^. Bound biexciton-based entanglement 
generation schemes [lO, Q, producing entangled polaritons with opposite spin, needs spe- 
cific tunings for efficient generation, and are expected to carry additional decoherence and 
noise due to scattering of biexcitons. Moreover linearly polarized single pump excitation (i.e. 
dp = cos 6*12;) -|- singly) = 2~"^/^(e*^|— ) + e^'*^|-|-))) cannot avoid the additional presence of the 
co-circular scattering channel which can lower polarization entanglement. The experimental 
set-up that we propose is a two-pump scheme under pulsed excitation, involving the lower 
polariton branch only. The pumps {pi and P2) are chosen with incidence angles below the 
magic angle so that single-pump parametric scattering is negligible. In this setup, mixed- 
pump processes (signal at in-plane wave vector k, idler at k, = ki + k2 — k) are allowed, as 
shown in Fig. [H As signal-idler pair, we choose to study the two energy-degenerate modes 
at k ^ {kix, k2y) and kj ^ {k2x, ^ly)- For all the numerical simulations we will consider the 
sample investigated in Ref. j9|. We choose ki = (0., 0.)/.im~^, and k2 = (0.9, 0.9)//Tn~^. In 
particular, we shall employ two pump beams linearly cross polarized (then the angle 9 will 
refer to the polarization of one of the two beams, see Fig. ([1])). This configuration is such 
that the counter-circular scattering channel (both bound biexciton and scattering states) is 
suppressed owing to destructive interference, while co-circular polarized signal-idler beams 
are generated. In the absence of the noisy environment, polariton pairs would be cast in the 
pure triplet entangled state |-|-, -|-) — exp (i4^)| — , — ). On the opposite side, two-pump exci- 
tation with counter-circularly pump beams (here not adopted) is able to stop this channel 
giving rise to the entangled state |-|-, — ) — exp (i46')| — , -|-). 

The advantages of this configuration are manyfold. First, detrimental processes for entan- 
glement as the excitation induced dephasing results to be largely suppressed [0. Spurious 
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coherent processes, e.g. Resonant Rayleigh Scattering [2l|, are well separated in k-space 
from the signal and idler modes. In addition, this configuration with signal and idler close 
;o the origin in k-space makes negligible the longitudinal-transverse splitting of polaritons 



22l | (relevant at quite high in-plane wave vectors). It is worth noticing that one could 
swap the roles of signal-idler and pump pairs, resulting in a parametric setup with isoener- 
getic pumps. However, such configuration would suffer from non-symmetric dephasing and 
non-equilibrated signal-idler photon intensities, that would lower the quantum correlation 
between the outgoing photons [23 1. 
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FIG. 1: (color online) (a) sketch of the proposed excitation geometry and of the lower polariton 
branch. The Gaussian pumps are linearly cross polarized with zero time delay. The specific 
polarization configuration with ^ = 0, where Q is the angle between e^i and the x-axis on the 
X — y plane, is depicted, (b) The simulated spectrally integrated polariton population in k-space. 
The parametric process builds up a circle passing through the two pumps and ^2)1 signal and 
idler polariton states are represented by any two points on the circle connected by a line passing 
through its center. For illustration the pair of signal-idler polariton modes chosen for entanglement 
detection are depicted as yellow crosses. The disc-shape contribution centered at the origin is the 
incoherent population background produced by phonon scattering. 



Following Ref. 



161 ] the tomographic reconstruction of the two-polariton density matrix 



is equivalent, in the o = {+, — } polarization basis, to the two-time coincidence 



(1) 



where Pl^^ (-Pk^a-) creates a signal polariton at k (an idler polariton at kj = ki -|- k2 — k) , A/" is 
a normalization constant and the detector window. We choose a very wide time window 
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Td = 120 ps, allowing feasible experiments with standard photodetectors. In order to model 
the density matrix eq. ([1]), we employ the dynamics controlled truncation scheme (DCTS), 
starting from the electron-hole Hamiltonian including two-body Coulomb interaction and 
radiation-matter coupling. In this approach nonlinear parametric processes within a third 
order optical response are microscopically calucated. The main environment channel is 



acoustic phonon interaction via deformation potential coupling 15|, |2J] . In order to account 
'or coherent and incoherent processes on an equal footing we use a DCTS-Langevin approach 
isl , with noise sources given by exciton-LA-phonon scattering and radiative decay (treated 

in the Born-Markov approximation). For mixed-pump processes with arbitrary polarization 

the Heisenberg-Langevin equations of motion read: 



In Eq. ([2]) "Pk o- are the projections onto the circular basis a of the coherent pump polariton 
fields. The complex polariton dispersion u)k includes the effects of relaxation and pump- 
ization (shifts), g is the nonlinear interaction term driving the mixed para- 



induced renorma! 

metric processes 15| ; summation over the repeated polarization indices ai and o"2 is assumed 
and the following selection rule holds: a + Uj = Ui + (J2. In general, third order contribu- 
tions due to Coulomb interaction between excitons account for terms beyond mean-field, 
including an effective reduction of the mean-field interaction and an excitation induced de- 
phasing. It has been shown that both effects depend on the sum of the frequencies of the 



scattered polariton pairs 



20|. For the frequency range here exploited, the excitation induced 



dephasing is vanishingly small and can be safely neglected on the lower polariton branch 



201 1 . On the contrary, the matrix elements of the (J = 0) counter-circular scattering chan- 



nel is lower (about 1/3) than that for the co-polarized scattering channel, but certainly not 



negligible 20|]. However, in the pump polarization scheme that we propose, performing the 
polarization sum in Eq. ([2]), it is easy to see that the counter-circular channel cancels out, 
as already pointed out. This feature is unique to the present scheme, while all previously 
adopted pump configurations suffer from the presence of both singlet and triplet channels. 
Equation ([2]) is a system of two coupled equations for polariton operators, acting onto the 
global system and environment state space, thanks to the two additive noise sources jFp^ , 
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FIG. 2: (color online) a(b) Real (Imaginary) part of the density matrix in the tomographic recon- 
struction according to eq. ([1]) in the linearly cross polarized pump configuration with 6 = 7r/6. 
This case results in E{p) ~ 0.7523 



Early experiments in semiconductor microcavities [13|, |27| provided promising, though 
indirect, indications of polariton entanglement. In order to achieve a conclusive evidence of 
entanglement one has to produce its quantitative analysis and characterization, i.e. a mea- 
sure of entanglement. Among the various measures proposed in the literature we shall use 
the entanglement of formation (EOF) E{p) 171] for which an explicit formula as a function 
of the density matrix exists. It has a direct operational meaning as the minimum amount 
of information needed to form the entangled state under investigation out of uncorrelated 



ones 
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The complete characterization of a quantum state requires the knowledge of its 



density matrix. For a quantum system composed of two two-level particles, it can be re- 



constructed by quantum state tomography |lol . [l6 |. It requires 16 two-photon coincidence 
measurements based on various polarization configurations [l^. Exploiting the Wick factor- 
ization 291] and the symmetries of the system, the density matrix elements are built up on 
the signal and idler occupation Ns/i± = {Pl^^^{T)Ps/i±{T)) and on the two-time correlation 
functions {P}_^{u)P-_^_{v)) and {P}_{u)P-_{v)). We integrate the system of equations eq. ([2]), 
coupled with the underlying non-equilibrium equations for the noise correlation functions 
(through time-dependent fluctuation-dissipation relations) 15|, l25|] . For a generic polariza- 



tion 9, the populations are independent on 6, N, 



s/i+ 



N. 



s/i- 



= Ns/ii whereas correlations 
satisfy {PIj^{u)PIj^{v)) = —e^^^{P}^{u)P-_{v)). In Fig. [2] a tomographic reconstruction is 
shown. We point out that different phase relations appearing in the non-diagonal terms of 
a reconstructed density matrix are directly related to the choice of pump linear polarization 
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FIG. 3: (color online) Dependence of the EOF on pumping intensity. The laser intensity / is 
measured in units of Iq = 21 photons /um~^/pulse according to Ref. Q]. 



6 (see caption of Fig. [2]). In the absence of longitudinal-transverse splitting 22j polariton 
entanglement is independent of the angle of the linear pump polarization 9. As Fig. [3] shows, 
there is a non-negligible region of the parameter space where, even in a realistic situation, 
high entanglement values are obtained. For increasing pump, EOF decays towards zero. 
This is a known consequence of the relative increase of signal and idler populations lOl. Ill] 
- dominating the diagonal elements of p - with respect to two-body correlations responsi- 
ble of the non-diagonal parts, which our microscopic calculation is able to reproduce. We 
expect entanglement to be unaffected by both intensity and phase fluctuations of the pump 
laser. The former are negligible according to Fig. 3, while the latter only act on the overall 
quantum phase of the signal-idler pair state. 

Different entanglement measures generally result in quantitatively different results for a 

B'ven mixed state. However, they all provide upper bounds for the distillable entanglement 
], i.e. the rate at which mixed states can be converted into the "gold standard" singlet 
state. Small EOF means that a heavily resource-demanding distillation process is needed for 
any practical purpose. Fig. [3] shows how a relative small change in the lattice temperature 
has a sizeable impact on entanglement. As an example, for the pump intensity / = 15/o, 
increasing the temperature from T = IK to T = 20/\ means to corrupt the state from 
E{p) ~ 0.88 to E{p) ^ 0.24, whose distillation is nearly three times more demanding. For 
a fixed pump intensity. Fig. [3] shows that, above a finite temperature threshold, the EOF 
vanishes independently of the pump intensity, i.e. the influence of the environment is so 
strong that quantum correlation cannot be kept anymore. In physical terms, at about 30K 
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the average phonon energy becomes comparable to the signal-pump and idler-pump energy 
differences, and the thermal production of signal-idler pairs is activated. 

Semiconductors heterostructures are complex systems in which other noise sources are 
expected. The simplest way to model this additional noise is via the introduction of a 
constant, temperature- and pump-independent, noise background A^;,. This quantity also 
accounts for the noise background characterizing the photodetection system. In the inset of 
Fig. [3], the dependence of the EOF on A''^ is highlighted. In our simulation, the quantity 
Ni, causes EOF to vanish in the limit of low pump intensity. As Fig. H] shows, a value of A^;, 



1 

P 



2 



of about 10""^ is realistic, as suggested by experiments [13| showing that it is considerably 
smaller than the PL-noise studied here. From inspection of the density matrix (in the 

(-|--|-, H — , — h, ) circular-polarization basis), for pump intensity / ^ 0, if the correlation 

dominates, we have as limiting triplet pure state 

/ 1 1 \ 


V 1 1 / 

on the other hand, if the population dominates, the state becomes separable p — > At 
leading order in J, the two contributions are comparable, and EOF emerges due to the 
normalization term A/". For nonzero A^;, instead, the leading term in the population is 
constant. This explain the behaviour of EOF as / ^ 0. Mathematically speaking, for any 
finite value of A";,, we have a separable state in the limit / 0. As seen in the inset of Fig. 
[3l a finite A";, affects only the range of very low pump intensity /. Figure 4 displays the 
time-resolved signal/idler occupations calculated at T = 5K for three different excitation 
densities. The dashed lines describe the PL contributions to the occupation. The figure 
clearly shows that at low excitation density the detected intensity in the signal/idler channels 
arises mainly from PL. Nevertheless the obtained EOF for these intensities is very high, 
contrarily to intuition, but in agreement with recent results This result puts forward 
the robustness of pair correlations and entanglement that can be evidenced, even when noise 
is the dominant contribution to one-body correlations. 

In conclusion we have shown that microcavity polaritons can be cast in an entangled state 
in a controlled way and we have given a ready-to-use realistic experimental configuration 
able to measure the EOF out of a dominant time-dependent noise background, without any 
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FIG. 4: (color online) Signal/idler occupations (dashed PL only, full PL+parametric) in time for 
three excitation densities {I = 5 Iq, I = 10 Iq and / = 15 Iq) for T = 5K are shown. The shape 
of the pump pulse is depicted for reference. 

need for post-processing. We point out that the outcome of a polariton parametric scattering 
is an entangled state of an electronic excitation of the semiconductor - the polariton pair - 
contrarily to parametric downconversion in a nonlinear crystal j^, where only the outgoing 
photons are entangled. Here, the emitted photons serve merely as a probe of the internal 
degree of entanglement. Thanks to their photon component, polaritons can sustain quantum 
correlations over mesoscopic distances inside the semiconductor. This is why they bear a 
unique potential as a controllable embedded mechanism to generate quantum information 
in a device and transfer it to localized qubits (e.g. spin qubits) over distances of microns j^. 

S.P. acknowledges the support of NCCR Quantum Photonics (NCCR QP), research in- 
strument of the Swiss National Science Foundation (SNSF). 
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